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Abstract: We show that brane inflation is very sensitive to tiny sharp features in extra 
dimensions, including those in the potential and in the warp factor. This can show up 
as observational signatures in the power spectrum and/or non-Gaussianities of the cosmic 
microwave background radiation (CMBR). One general example of such sharp features is 
a succession of small steps in a warped throat, caused by Seiberg duality cascade using 
gauge/gravity duality. We study the cosmological observational consequences of these steps 
in brane inflation. Since the steps come in a series, the prediction of other steps and their 
properties can be tested by future data and analysis. It is also possible that the steps are 
too close to be resolved in the power spectrum, in which case they may show up only in 
the non-Gaussianity of the CMB temperature fluctuations and/or EE polarization. We 
study two cases. In the slow-roll scenario where steps appear in the inflaton potential, the 
sensitivity of brane inflation to the height and width of the steps is increased by several 
orders of magnitude comparing to that in previously studied large field models. In the 
IR DBI scenario where steps appear in the warp factor, we find that the glitches in the 
power spectrum caused by these sharp features are generally small or even unobservable, 
but associated distinctive non-Gaussianity can be large. Together with its large negative 
running of the power spectrum index, this scenario clearly illustrates how rich and different 
a brane inflationary scenario can be when compared to generic slow-roll inflation. Such 
distinctive stringy features may provide a powerful probe of superstring theory. 
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1. Introduction 



It is generally believed that the early universe went through an inflationary epoch In 
the field theory context, the most studied realization is the slow-roll inflationary scenario 
Q. Such models confront primordial cosmological observation via two data points : the 
magnitude of the density perturbation (or temperature fluctuation in the cosmic microwave 
background radiation (CMBR)) as measured by COBE ~ 10~^, and the approxi- 

mately scale invariant scalar power spectrum, with index n<j ~ 1. In addition, there are 
three observational upper bounds on, namely, the tensor to scalar ratio r < 0.3, the run- 
ning of the spectral index \dns/dlnk\ < 0.1 and the non- Gaussianity /nl < 0(100) @,§. 
Typical predictions of dus/dlnk and /nl (and maybe even r) may very well be far too 
small to be detected. Since the COBE normalization typically fixes the scale of the model, 
so far, the deviation of from unity is the only parameter that discriminates among slow- 
roll inflationary models. It is not surprising, therefore, that present data allows a large 
degeneracy of models . 

If superstring theory is correct, it must provide an inflationary scenario for the early 
universe. It should provide a precise and restrictive framework for inflationary model 
building while at the same time, new scenarios and new phenomena are likely to appear. 
Hopefully it will predict distinctive stringy signatures that, if observed, can be used to 
support superstring theory. It is along this direction that this paper attempts to address. 
Here we like to discuss the possibility of a stringy feature that may show up in a distinctive 
fashion in the CMBR. If the proposal is realized in nature, the number of data points 
linking theory and cosmological observation may increase (from ~ 2) by as much as an 
order of magnitude, substantially strengthening the test of inflation as well as the specific 
stringy realization of inflation. 

It has become clear that the brane inflationary scenario is quite natural in string theory 
0. Here, the inflaton is simply the position of a mobile brane in the compactified bulk 
||8|, In flux compactification in Type IIB string theory, where all closed string moduli 
are dynamically stabilized, warped geometry appears naturally in the compactified bulk 
p!o| ]. Warped geometry in the 6-dimensional compactified bulk tends to appear as warped 
throats. In the simple but realistic scenarios, inflation takes place as D3-branes move up or 
down a warped throat |11, 12, IS]. Inflation typically ends when the D3-branes reach the 
bottom of a throat and annihilate the anti— D3-branes sitting there. The energy released 
heats up the universe and starts the hot big bang. 

The warped geometry either shifts masses to the infrared and helps to flatten the 
inflaton potential, or provides speed-limit that restricts the inflaton velocity, thus improving 
the e-folds of inflation. A typical throat is a warped deformed (or resolved) conifold. A 
string scale test mass in the bulk is red-shifted to a smaller value at the bottom of the 
throat. Via gauge/gravity duality, moving down the throat corresponds to a gauge theory 
flowing towards the infrared region. The best known example suitable for brane inflation 
is the M = 1 supersymmetric Klebanov-Strassler (KS) throat in type IIB string theory 
on approximate AdS^ x T^'^ background |l^. It is argued that the KS throat is dual to 
the TV = 1 supersymmetric SU{N + M) x SU{N) (with N = KM) gauge theory with 
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bifundamental chiral superfields and a quartic tree level superpotential in four dimensions 
m, [l^. As the gauge theory flows towards the infrared, SU{N + M) = SU{{K + 1)M) 
becomes strongly coupled and Seiberg "electric-magnetic" duality says that it is equivalent 
to a weakly coupled SU{{K — 1)M) with appropriate matter content |17]. As a result, the 
gauge theory becomes SU{{K — 1)M) x SU{KM) with the same set of bifundamentals. 
This process repeats as the gauge theory continues its flow towards the infrared; that is, it 
goes through a series of Seiberg duality transitions. This is known as the duality cascade 

Hi- 

Recently, a more detailed analysis of this Seiberg duality cascade [^] shows that there 
are corrections to the anomalous mass dimensions and their effects in the renormalization 
group flow of the couplings on the gauge theory side show up as small structures on the 
supergravity side: as steps appearing in the warped metric. The steps in the warp factor are 
expected to be slightly smoothed out, leading to a cascade feature that may be observable 
in cosmology. It was conjectured in Ref.|19] that the origin of the steps might be localized 
charges (branes) at radial positions inside the throat. To fully understand the gravity dual 
of the duality cascade, one needs to take the branes into account and study the resulting 
supergravity solutions together with the thickness and the stability of the branes including 
other smooth corrections to the KS solution. We leave this issue for future study and in 
this paper we focus only on the observable effects due to the steps. The presence of steps 
is quite generic: either as a correction to the approximate geometry for a warped throat, 
or as a variation of the simplest warped deformed throat in the KKLMMT-like scenario. 
Although the KS throat is the best understood example, in some sense it is too simple, i.e., 
there are no chiral fields in the infrared. Warped throats are natural consequences of flux 
compactification. They arise by considering Calabi-Yau singularities that admit complex 
deformations, corresponding to smoothing the singular points with 3 cycles on which to 
turn on fluxes. The holographic duals are given by duality cascades of gauge theories on 
D3-branes at the singular geometries. Generically, steps due to jumps in magnitudes of 
anomalous mass dimensions appear. Each duality cycle contains a number of steps (the 
KS throat has only one step in each cycle). Such steps can introduce sharp features in 
the cosmic microwave background radiation. Such a step may have been observed already 
1^, P]. In general, a D3 brane encounters many steps as it moves inside a throat. If 
the steps are close enough, as expected, the power spectrum may be changed little, but 
the impact on the non-Gaussianity can be large. The duality cascade behavior of the 
warped geometry can therefore have distinct observable stringy signatures in the cosmic 
microwave background radiation. The possibility of detecting and measuring the duality 
cascade is a strong enough motivation to study the throat more carefully. Although both 
Seiberg duality and gauge/gravity duality are strongly believed to be true, neither has been 
mathematically proven; so a cosmological test is highly desirable. This will also provide 
strong evidence for string theory. 

It was proposed that one can easily envision a multi-warped throat, that is, a throat 
inside another throat (i.e., the IR of one throat is matched to the UV of another throat) 
|2C]. Here we expect a kink in the warp factor, but not a step. Again, we expect this to 
generate some non-Gaussianity. 
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Here we discuss the generic features introduced by the presence of such steps in the 
warped geometry of a throat. We consider the presence of steps in both slow-roll and 
DBI inflation. To be general, we shall allow the steps to go either up or down and we 
consider both cases. Their properties can be quite different in these scenarios. However, 
some features can be quite independent of details. For example, as an illustration, let us 
assume the feature in WMAP CMBR data at / ~ 20 (with width about 61 ~ 10) is due to 
such a step. If we further assume that the deviation around / ~ 2 is at least partly due to 
another step (i.e., not all due to cosmic variance), then we can predict that its width SI ~ 2. 
Furthermore, we would then expect additional steps at / ~ 200 (with width 61 ~ 100) and 
/ ~ 2000 (with width 61 ~ 1000) and so on. However, the sizes of the steps are more model 
dependent, though they tend to vary monotonically. If there is no step at / ~ 2, then the 
next step is expected at some / > 200. 

As illustrations, we consider two specific scenarios: 

1. Its properties in slow-roll brane inflation. Effects of steps on power spectrum in 



CMBR in slow-roll inflation have already been studied in inflationary models [21, 22 



24] and in details |25, 26, 27| in chaotic inflation. Its effect on the non-Gaussianity 



has also been studied in details in Ref. |28, Assuming that the / ~ 20 feature 



is due to a step, Ref. ||2^ shows that the non-Gaussianity in chaotic-like inflation 
gives 

^feature ^ o{lO). Because the shape and running of this non-Gaussianity is 
distinctively different from the cases done in data analysis, whether the PLANCK 
satellite is going to be able to test it remains an open question. In slow-roll brane 
inflation, the inflaton, being the position of the brane, is bounded by the size of 
the flux compactification volume, and so typically (p <C Mpi, as opposed to large field 
case for chaotic inflation, where the inflaton takes values much larger than the Planck 
mass {(p > 15Mpi). Consequently, the sensitivity of brane inflation to the height and 
width of the steps is increased by several orders of magnitude. 

However, the appearance of the warp factor in the potential typically destroys the 
flatness of the potential, so some fine-tuning is required. Here, we shall simply assume 
slow-roll with steps in the small field case. 

2. Properties of steps in the infrared Dirac-Born-Infeld (IR DBI) model. In the UV 



DBI model, Z)3-branes move down the throat at relativistic speed |12, In terms 



of known brane inflation scenarios, this model produces too large a non-Gaussianity 



to satisfy the current observational bound [^, |3^, |3J] and inconsistent probe 
brane backreactions |^5|. In the IR DBI model, inflation happens when L)3-branes 
are moving relativistically out of a throat |13, 36]. Recent analysis ]^5[ shows that 



the IR DBI brane inflationary scenario fits the CMBR data quite well, comparable 
to a good slow-roll inflationary model. However, it has a regional large running of 
n^: —0.046 < dus/dlnk < —0.01 at 0.02/Mpc due to a stringy phase transition 
in primordial fluctuations. The spectral index has a red tilt ~ 0.94 at /c ~ 
0.02/Mpc, and transitions to blue at larger scales. It also predicts a large distinct 
non-Gaussianity due to the DBI effect, which should be observable by PLANCK. So 
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we see that such IR DBI models of brane inflation can flt the present cosmological data 
and have distinctive stringy predictions difi'erent from slow-roU. In this DBI scenario, 
steps appear in the warp factor. The feature appearing in the power spectrum has 
different properties from that in the slow-roll case. For example, we can choose the 
phenomenological parameters that give reasonable fit to the glitch around / ~ 20, but 
given the extra parameters that are added, the improvement in is not significant. 
We will explain why, in general, the glitches caused by the warp factor steps in IR DBI 
model is much smaller than the glitches in slow-roll inflation, while the associated 
non-Gaussianity signal can be very large. In this model, the non-Gaussianity should 
come from a combination of DBI relativistic roll effect together with the step effect, 
which can be large and is very distinctive. 

As we shall see, some of the predictions are rather sensitive to the details, while others 
are robust, for example, that there are a number of steps, with their spacing A In oc oc 
A In / roughly equal. 

In the IR DBI model with steps, CMBR data can provide the following data points: n^, 
dug/dln k, and DBI non-Gaussianity, which may be crudely quantifled into four data points 
(the magnitude as measured by Jnl, two data points from the shape of the bi-spectrum, 
and another for its running^). So the IR DBI scenario alone may confront CMBR via six 
data points. Comparing to the slow-roll case, which so far confronts CMBR via only n^, 
the test of IR DBI is clearly much more restrictive. Note that the number of microscopic 
parameters of the model is roughly flve (depending on how one counts) in both brane 
inflationary scenarios. Now, let us include the steps. Each step will provide three data 
points: its position, height and width, in addition to the non-Gaussianity associated with 
the step. The number of microscopic parameters of the model will increase by at least 
three. If we see more than one step in the CMBR, the test of the underlying superstring 
theory will be very stringent. In fact, we can learn a lot about the flux compactification 
of our universe in superstring theory. Of course, one may add the cosmic string prediction 
into the list to further test the scenarios. 

We should mention that, broadly speaking, there are two classes of inflationary sce- 
narios in string theory, depending on whether the inflaton is an open string mode (brane 
inflation) or a closed string mode (moduli inflation)^. Moduli inflation can reproduce usual 
slow-roll inflation, but so far no distinct stringy signature has been identifled. Obviously it 
is important to continue the search for distinct stringy signatures in both brane inflation 
and moduli inflation. 

This paper is organized as follows. Section ^ reviews the setup of the geometry and 
the model. It also reviews DS-brane potential where the step appears as a small correction. 
Section ^ presents general formalism on the power spectrum and bispectrum due to a sharp 
step. Section Q discusses the step in the slow-roll brane inflation. Section ^ considers a 

^This accounting is clearly an oversimplification, since the bi-spectrum is a function. Furthermore, the 
4 point-correlation (tri-spectrum) has also been calculated which may be measured and tested. 

^In brane inflation, infiatons are positions of branes, which are vacuum expectations of light open string 
modes. However, the infiaton potential can still be described by closed string interactions, perturbative, 
collective and/or otherwise. 
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step in the IR DBI model. We perform qualitative analyses of the step features in IR DBI 
inflation followed by numerical analyses. The comparison to WMAP data around I ~ 20 
is also performed. To be more general, Sec. ^, ^ and |5| are written in a way that the 
main analyses apply to general sharp features, such as steps. The duality cascade steps 
reviewed in Sec. ^ is used as an example to give a more restrictive set of predictions. As 
we emphasized, some details of the steps such as the width, depth and forms may vary 
from example to example. In Section |6| we make some remarks on some broader aspects of 
brane inflation. A number of details are contained in the appendices. 



2. The Model 



First it is useful to define a basis of coordinates and a metric. The ten dimensional metric 
which describes the throat is that of AdS^ x X^, where has the T^'^ geometry in the 
UV region. Including the expansion of the universe, the metric has the form 



h\r){-dt^ + a(t)2dx2) + h'\r){dr^ + r^ds'j^^), 



(2.1) 



(is^i,i is the metric that describes a base 



Far away from the bottom of the throat, dsx; 
of the conifold T^'^ which is an S"^ fibered over 5^. Here h{r) is the warp factor, that is, a 
generic mass m mh{r) in the presence of h{r). 

The inflaton (p is related to the position of 4-dimensional space-time filling D3- 
branes. In IR DBI inflation, they are moving out of the B throat into the bulk and then 
falling into the A throat. Inflation takes place while the branes are moving out of the B 
throat. In UV DBI inflation or the KKLMMT scenario, inflation takes place when they are 
moving down the A throat. Specifically, we use a single coordinate to describe collectively 
the motion of the branes, 

</> = ^/nsT-ir . (2.2) 
The edge of the throat is identified with = ^/n^T^R, with 



27 



i?* = —T^g^Na''^, N = KM 



(2.3) 



In regions where several kinds of backreactions from the 4-dimensional expanding back- 



ground [36, 38|, the stringy effects |31, pq] and the probe-branes [12, 39[ can be ignored 



the following DBI-CS action describes the radial motion of the D3 branes. 



S 



I 



d x^P^ 



/. 5"^^^'^^^'^+r(,^)_y(^) 



-e-*r(</>)Ji 



The warped Z)3 brane tension T((/)) and the inflaton potential is given by 

r(0) = UBTsh^ct)) , 



V, 



DD\ 



Vn 1 



47r2v (0 - (j)Ay 



(2.4) 

(2.5) 
(2.6) 

(2.7) 
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Here v is the volume ratio; one may take v = 16/27. For our purpose here, the Coulomb 
term is negligible so V^fj{(j)) = Vq = 2nAT3h'^{(j)A) is the effective vacuum energy. We 
expect 1/3 1 ~ 1. Positive /? » 1 for UV DBI model while negative /3 ~ — 1 for IR DBI 
model. Slow-roll inflation requires |/?| <C 1. 

It is convenient to introduce the inflationary parameter e = —H/H'^, so that 

- = H\l-e). (2.8) 
a 

The universe is inflating when e < 1. In all the scenarios that we consider here, e never 
grows to 1, so inflation ends by DD annihilation. 

For simplicity, let us assume that the B throat is a KS throat. The warped factor for 
the B throat around the pf^ duality transition (starting from the bottom of the throat) is 
simplified to (see Appendix ^ 



h (r) ^ — (1 + A) , A = ) 



R%Pi ' ^ 16vr p3 



1 + tanh 



dp 



(2.9) 



where h{r = Rb) — 1 at the edge of the throat. Here rp is the positions of the steps, the 
initial 3> 1 so that the warped factor formula is approximately good, and dp controls 
the width of the step. The steps has a constant separation in Inr, that is, Inrp+i — Inr^ ~ 
2Tr/3gsM. As one moves down the throat (r and p decreasing), note that the step in the 
warped factor h^{r) of the KS throat is going down. See Figure |l[ This stepping down 
happens at each Seiberg duality transition. Together, they form a cascade |^8|. For large 
K, one encounters K — 1 steps as one approaches the infrared. Note that we are ignoring 
smooth corrections to the shape of the warp factor even though they may be larger than 
the step size. This is a reasonable approximation since it is the sharp features that will 
show up as distinctive features in the CMBR. 

Steps are a generic feature: consider another throat whose gauge dual has uq gauge 



factors, with appropriate bi- fundamentals |2C]. As the gauge coupling of the gauge factor 



with the fastest running towards large coupling in the infrared (decreasing r) gets strong, 
Seiberg duality transition applies. This happens for each gauge group factor sequentially 
until we reach the same structure as the original gauge model. Typically, it takes nc/n^ 
number of transitions to complete this cycle, where Us is a factor in nc (the KS throat has 
no = n-s = 2 so each cycle has only 1 transition). Because of the jump in the corresponding 
anomalous mass dimension, steps in the warp factor are very generic. It is also likely that 
for another throat with a different geometry, some steps may show up with an opposite 
sign. To be general, we shall discuss each of these scenarios. If the throat has relatively 
few steps, or if the steps are well separated, the impact of individual step on CMBR may 
be observable. Otherwise, the steps may be too close for them to show up in the power 
spectrum. 

In the original KS throat solution, they take the approximation where the dilaton is 
constant. Here the dilaton factor, e~*, runs in general and is (j) dependent. However, for 
the IR DBI inflation model we discuss in this paper, most of the DBI e-folds are generated 
at the tip of the throat where HR'^/Nf, < r < HR^. At the same time, the brane moves 
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Figure 1: The warp factor h{r) in the KS throat, including the steps. Here r is measured in \/a' 
with Rb =i 100, gs = 2, M=20,K = 10 and N = 200. The width d = IQ-^. In this figure, there 
are actuahy 4 steps, located at r ~ 73, r ~ 82, r ~ 89 and r ~ 95, although the step at r ~ 95 
is too small to show up. Here, the parameters (in particular, a large gsM) are chosen so that at 
least 3 steps are big enough to show up in the figure. This leads to relatively large corrections 
to the positions of the steps. Other parameters should be used in more realistic situations and in 
comparison with data. 



across roughly gsM steps [35|. The running of the dilaton can be ignored if ggM <C pi, 
which is most easily satisfied with a small gs- Furthermore, the WMAP data only covers 
a few e-folds, which corresponds to 3gsM/Ne steps, so we only need pi > ggM to safely 
ignore the dilaton modification to the kinetic term. In this paper, when we discuss the IR 
DBI inflation, we always absorb into gs as a constant, so it will never appear explicitly 
in our analysis. 

The energy density and pressure from the DBI action are given by 

p = V{4>) + T{(t>){ci' -1), p = -V{(t>) + T{<p){l-Cs), (2.10) 
where the sound speed is defined as the inverse of the Lorentz factor 7, 



Cs = 7"' = V 1 - '^'/r • (2.11) 
The background equation of motion is given by 

V{(t>)+T{4>){c;' -l) = 3H\ (2.12) 

^ - l^'}^' + ^HcU + T'i<l>) + cliV'i^) - r'(0)] = . (2.13) 

On the other hand, the running dilaton is more relevant for the UV DBI or the 
KKLMMT scenario, because here infiationary e-folds are generated by the D3 brane mov- 
ing all the way from the UV to the IR end of the throat. For the KKLMMT scenario, we 
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expand the DBI action in the cj)^ <^ T{cl)) hmit, and get the slow-roll action, 

^1 



SsR = d X 



e-*<A2-r(</.)(e-*-i)-y(0) 



(2.14) 



We see explicitly that the effective inflaton potential, VcQ{<j)) = r(0)(e~* — 1) + V{(t)), 
now exhibits the step feature from T{(l)). In the original KS solution, e~* = 1, and the 
steps in the warp factor do not affect the inflaton potential. The term T((/))(e~* — 1) 
is generically bad for slow-roll inflation, however. Since T((/>) ~ (/>^, we have effectively 
introduced a quartic term in the potential. It is well known that a Ac^^-type potential 
requires a trans-Planckian (j) to generate inflation. However, in the KS throat (and other 
GKP type compactifications |]lO[), trans-Planckian (p is impossible ^ and we have a 
strong bound on the largest inflaton field value, 

In order to get slow-roll inflation, the dominant term in the potential has to be the uplifting 
term from D3 at the bottom of the throat, so we require 

yo»T(</.)(e-*-l) . (2.16) 

Since Vq ~ h'^{(/)A) *C T{(p) ~ h'^{(p), the above relation is true only when e~* ~ 1. Note 
that if Vq dominates the inflaton potential, the step height in the inflaton potential is not 
given by A as defined in Eq.( |2.9D , but is suppressed by a factor (e~* — l)(/)'^/4>\ <C 1, i.e. 

^ = (e--l)^A. (2.17) 

So far, we have argued that even if the dilaton runs generically, slow-roll inflation can only 
occur within the region where e~* stays close to 1. Thus when we discuss slow-roll scenario 
in this paper, we will ignore the e~* modification to the kinetic term, and set e~*0^ « <j)'^, 
but we keep the term T(0)(e~* — 1) in the potential. This approximation should be pretty 
good, since e~* only exhibits mild features like kinks, while T{(l)) has sharp features like 
steps. When both of them are present, the sharp features in T((j)) definitely dominate. 

3. The Power Spectrum and Bispectrum in General 

The gauge invariant scalar perturbation (^(r, k) can be decomposed into mode functions 

C(r, k) = u{t, k)a(k) + u*{t, -k)a^(-k) . (3.1) 
The mode function Uk(''") is given by the equation of motion 

v'l + (^k'^ci -^^vk = 0, (3.2) 

where 

Vk = zuk, z = a\f2elcs, (3.3) 
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and prime denotes derivative with respect to the conformal time r (defined as dt = adr). 
In the absence of any features in either potential or warp factor, the power spectrum is 
given by 

Pn{k) ^ l^kkP (3.4) 



(3.5) 



where Uk is u{t, k) evaluated after each mode crosses the horizon. 

To study the effects of features, we define the inflationary parameters as below 

^=-^' ^=^' ^=^- 

We can express z" / z exactly as, 

--2aH + + _ + + _ + . (3.7) 

The z"/z encodes all the information from the inflationary background, and determines 
the evolution of ti(T, k). In the absence of sharp features, e, ry and s remains much smaller 
than 1, so z" /z ~ 2a'H'. However, a sharp feature in the inflation potential or the warp 
factor will induce a sharp local change in e, fj and s, and z" /z has a nontrivial behavior 
deviating strongly from 2a?H^ around the feature. In this paper, we have analyzed two 
cases. First, the slow roll brane inflation with a step feature in the inflaton potential V{(j)). 
Due to the small field nature of brane inflation, e is negligible, which greatly simplifies z" / z 
to (See Appendix |B|) 

£:.....(._]^).. 3^). 

Second, the IR DBI brane inflation scenario with a sharp step appearing in the warp 



factor. In Section 5.2.2, we show that the sharp change of the sound speed is the major 



contribution to z" / z, and we have 



where T is the warped brane tension defined in Eq.(p.5|). To avoid confusion, we emphasize 
that the primes on z or Uk (functions of time) denote the derivative with respect to the 
conformal time r, while the primes on T or y (functions of (j)) denote the derivatives with 
respect to the field (j). 

Among all the terms in the exact cubic action for the general single field inflation |pT|l , 
the 3pt for sharp features receives dominant contribution from the term proportional to 
(t-) in most interesting cases. 



dr 

>3 



(27r)3<53(^k,) + c.c. . (3.10) 
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The reason is a generalization of that in Ref. |28]. As we wih show, e and Cg do not change 



much due to approximate energy conservation relations; but they change in a very short 
period of time, boosting e or Cg by several orders of magnitude. 

4. Steps in Slow-Roll Brane Inflation 

In the slow-roll scenario, let us assume that inflation takes place as the brane moves down 
the throat. Since the dilaton runs, the slow-roll potential is given by 

Vsr{c^) = V{ct^) + r(0)(e-*W - 1) . (4.1) 

Here, the second term in the potential does not vanish, so the step in T{(j)) shows up in 
the potential. Again, we may consider both step up and step down cases. 

Before we proceed, we mention two other possibilities that will not be studied in more 
detail in this paper. Suppose the warp factor steps down as r increases. As the inflaton 
moves down the throat, it encounters such steps. If the height of any particular step is 
too big, then the inflaton cannot overcome it so its classical motion towards the bottom 
of the throat is stopped. Let us say this happens at r = Vp. It can still move along the 
X S"-^ angular directions. In a realistic model, the x 5^ symmetry is slightly broken 
so the inflaton would move along the steepest angular direction, until it reaches the lowest 
point at this radius rp. If the step there is lessened so that the inflaton can roll over it, the 
inflaton can continue to roll down along r until it reaches the next step. This process can 
repeat some number of times. As a result, the inflaton path is substantially longer than 
just moving down along r. Effectively, this may substantially extend the field range and 
enough e-folds of inflation is much more likely. 

Another possibility is that the brane can tunnel over the step barrier. Recent analysis 
of the tunneling of branes suggests that, under the right condition, the brane may tunnel 



very easily across a potential barrier that is not particularly small |42]. If this is the case, 
then the brane may roll along some angular direction until it reaches such a position and 
tunnel right over the potential barrier there. It can then roll down along r until it reaches 
the next step. Again, this process can repeat any number of times, as a result, the inflaton 
path is again substantially longer than just moving down along r, allowing the generation 
of many e-folds. 

In the following, we consider the scenarios where the inflaton can roll over many steps 
in the potential. 

4.1 The Position of the Steps 

The duality cascade predicts more than one step, so it is interesting to ask where the steps 
are. Assume there is a step at the / ~ 20 in the CMBR, which may be checked with a non- 
Gaussianity measurement. Let us further assume that there is another step at / ~ 2, (the 
significance of the signal at / = 2 in power spectrum is much smaller than that in / ~ 20 
due to the large cosmic variance, so we assume this mainly for the purpose of illustration). 
We now show that the location of all the steps in the duality cascade can be determined 
by given the locations of two such adjacent steps. 
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Using I ~ 10^(A:/Mpc~i), we find 

- diVe din ~ din/ ~ Flit ~ , (4.2) 

Since both H and (p are slowly varying during inflation, we have 

dlnl(xd4>. (4.3) 

Suppose 4> is decreasing (going down a throat), also suppose that the step at / = 2 is at 
(prno s-ncl that at / = 20 is at (pmo+i, we have 

0mo+l 4'mo+2 

For large ^^M, this ratio is close to unity, e'^'^^^^"^ ~ 1 + so <pmo ~ 4>mo+i — 'Pmo+i ~ 
4'mo+2 — <Amo+i<^- Due to ([4.3D, equal spacing in cp implies equal spacing in In/. So we find 
that the next two steps are at around / ~ 200 and / ~ 2000 respectively. In addition, the 
effect of the step at multiple moment / should span over Al multiple moments with A/ oc /. 

4.2 The Power Spectrum 

In the slow-roll scenario without features, we have the attractor solution, in units of Mp, 

= V{cP)/3, 3H4> = -V'{cf>) , (4.5) 

where / is derivative with respect to the inflaton (j){t) and dot is derivative with respect to 
time. Here, 

H /H'V 1 ^V' 



It is convenient to introduce another inflationary parameter 

f, = e/He = -2r]sR + 4€sr (4.7) 
where the usual slow-roll version is given by 

V" 

rjSR = -y. (4.8) 

This yields n<j — 1 = — — 2e = 2rjsR — QesR- We emphasize that the above relations 
between e, fj and esR, rjsR, only hold for the attractor solution in absence of sharp features. 
We will use the parameters e and fj in our analyses for the sharp feature case. 

To see the full details of the effect of a step in potential, numerical calculation is 
necessary. However, the qualitative behavior can be estimated as follows. The step in 
the potential is typically characterized by two numbers: the depth which we describe as 
the ratio AV/V ~ 2c, and the width A0 = 2d in unit of Planck mass. We can divide the 
motion of inflaton into two parts: acceleration and relaxation. First, the inflaton, originally 
moving in its attractor solution, momentarily gets accelerated by the step. The potential 
energy AV gets converted to kinetic energy. Then after the inflaton moves across the step. 
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it starts to relax back to its attractor solution under (j) + 2>H(j) ~ 0, where the Hubble 
friction term dominates over the potential. 

We first look at the power spectrum. The inflaton velocity in the original attractor 
solution is given by \4>attr\ ~ |F'|/3-ff \/eVj?>. After the acceleration of the step, (j) ~ 
^V{2>c + e/3). (It turns out that the best fit data give comparable 3c and e/3, so for the 
purpose of order-of-magnitude estimate we later will also approximate (j) to be ^/2>cV .) A 
useful formula to understand the dip of the glitch in the power spectrum caused by the 
step is P-ji = H"^ / {Att"^ (fP') . The acceleration of inflaton does not change the H so it reduces 
the P-ji- The ratio of Pti between the initial value and the value at the first dip is related 
to the corresponding ratio of (j), which is y^l + 9c/e. For the large field quadratic inflation, 
Ref. ^ give the best-fit data, c = 0.0018, e = 0.009. As we will see later, the values of 
c and e will be rather model-dependent. However the discussion in this paragraph shows 
that their ratio should be fixed for a specific observed feature, 

- f« 0.2 . (4.9) 

e 




Figure 2: In the small field case, we show how the step in the inflaton potential changes the power 
spectrum. We show the power spectrum in the left panel and the behavior of e around the step in 
the right panel. For illustration, we use the KKLMMT scenario with only the Coulomb potential. 
We choose the background flux TV = 2000, so A</)/Mpi < 0.01. The inflation scale Vq ~ lO-^Afj^j. e 
is tiny, typically e ~ 10~^^ as shown in the right panel. We introduce a step with c = 8 x 10^^^ and 
calculate the power spectrum for four different values of d: (1) d = 3 x lO^^Afpi, green solid line. (2) 
d = 1.7 X lO-'^Mpi, black solid line. (3) d = 1.7 x lO-'^Mpi, red dashed line. (4) d = 1.7 x lO^^Mpi, 
blue dotted line. Note that the dip in the power spectrum depends weakly on d. In case (1) and (2), 
Ae has not saturated the bound Eq. 4.1C , so decreasing d will enhance the bump in Pn significantly. 
In case (3) and (4), where Ae is maximized, the bump in Pqi, does not depend sensitively on d, 
but the range of the oscillations in fc-space does. The black solid line is close to the best fit power 
spectrum given in Ref. ||2^, |2^ . 
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Although the above discussion indicates that the depth of the dip is relatively indepen- 
dent of d, the peak of the first bump is. When the inflaton moves across the step, energy 
conservation requires that ^A(/)^ < /SV , with the bound saturated when the Hubble friction 
is negligible. This will result in an upper bound in the change of e 

Ae w /W/H"^ < 5c . (4.10) 

We expect that the effect of decreasing d will first enhance the bump, since smaller d leads 
to more deviation from the attractor (larger Ae), and will enhance the bumps in Pji during 
the relaxation period. However, as long as Ae saturates the upper bound (when d is small 
enough to ignore the Hubble friction), further decreasing d will not affect the bumps in 
P-ji any more, since the relaxation period essentially starts with the same <j) no matter how 
small d is. So we expect the bump to depend sensitively on d if moving across the step 
takes 0(1) e-folds, and it becomes relatively insensitive to d when the step is so sharp 
that moving across it takes only <C 1 e-folds. In the latter situation, reducing d will only 
increase the extension of the oscillations in Pr, 

Akr^Mr^ko^, (4.11) 



where ko is the starting point of the feature, and (|3.8| ) is used. 



Ref. 1 26, 27] fit the / ~ 20 feature in the WMAP data introducing a step feature in the 
large field slow-roll model and gives the best fit power spectrum. Here we have numerically 
reproduced a similar power spectrum for the small field case. We are not performing a 
complete data analysis to find the best fit model here, our major emphasis is to show how 
the width d comes into play. Fig. |^ shows the power spectrum and e for four different 
values of d. The black solid line represents the power spectrum close to the best fit model 
given in Ref. (2|, with c = 8 x 10"^^ we find that d = 1.7 x lO^^Mpi. 



4.3 Non-Gaussianities 

Now let us look at the bispectrum. For slow-roll inflation, the 3-point function of the scalar 
perturbation (^(rend, k) includes terms proportional to e^, e^, and e-q' [43|. For the slow-roll 



potential without any features, the first terms dominate. Comparing to WMAP's ansatz, 
it gives the non-Gaussianities estimator /tvl = 0(e)- 

In the presence of sharp features, from ( [4.10| ) we see that e still remains small, < 
0(0.01). But fj' can be much larger. Therefore, the leading three-point correlation function 
is given by the term proportional to ef]' [p^ j, 

{C{Tend, kl)C(T-end, ^2)CiTend, ka)) 

= i ^Yluk,{Tend)^ J dra^ ef]' (^uI^{t)uI^{t)-^uI^{t) + peim^ 

X {2TTf6^{Y,^i)+c.c. , (4.12) 

i 

where the "perm" stands for two other terms that are symmetric under permutations of the 
indices 1, 2 and 3. The details of such an integration are quite complicated. Nonetheless we 
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can estimate the order of magnitude of the non-Gaussianity estimator fi\iL by comparing 
it to the slow-roll case. The most important difference is that here we replace a factor of 
e by r/', and also fj' only gets large momentarily. Hence we estimate f^f^^'^'^ = 0{f]'AT) = 
0{Afi), where At is the conformal time that the inflaton spends crossing the step. 

To estimate the level of these non-Gaussianities, we now give a qualitative estimate 
for fj' |28]. During the acceleration period, the e is increased by 



Ae w AV/H^ w 5c . (4.13) 

The duration of this period is 

Ataccel ~ A(^/0 « d/VcV , (4.14) 



where we used d) estimated in Sec. 4.2. These can be further used to estimate 



Ar? « = — « , (4.15) 

The time scale for the relaxation period is of order H~^, during which fj is 0(1) and ry are 
of order 0{H). 

We sum over the contributions from both the acceleration and relaxation periods. The 
former gives ffff^^ ~ 7c^/^/((ie), the latter gives = So for most interesting 

cases where non-Gaussianities are large enough to be observed, it can be estimated by the 
first contribution, 

7„3/2 

/feature ^ __ (4_ig) 

We note that this is only a crude order of magnitude estimation on the amplitude of this 
non-Gaussianity, since the integration ( 4.12| ) also involve mode functions which will 



be modulated by the presence of the sharp feature, and the shape and running of such 
non-Gaussianities are very important. Details have to be done numerically, as in 



Qualitatively since we have argued in Sec. 4.2 that for a specific observed feature the c/e 



is hold fixed model-independently, ( [4.16| ) implies that the value of d is very crucial to the 
level of the non-Gaussianities. There are two major constraints on d. First, in the power 
spectrum, as shown in Fig. § the bump in Pn may depend sensitively on d. Second, the 
range of oscillation in PTi{k) is also controlled by d. If d is too small, the oscillation in 
PTz{k) might spread over to the well measured first acoustic peak in WMAP Q curve. 



Numerically, we have found that y/c/d ~ 0(1), consistent with ( 4.11 ), so the magnitude of 
( |4T6|) should be close to that in Ref. |2|]. 



It is instructive to split this expression to fj^'^'^^^ ~ 7c/e • ^fcjd. The first factor 
is determined by the amplitude of the glitch from ( |4.9| ), while the second factor is the 
extension of the glitch from ( 4.11| ). Note both are in the /c-space not the CMB multipole 



/-space. Therefore in principle a sharp feature can also appear only in the 3-point function. 
This is clear from our estimation (|4.9| ) and ( [4.16| ), where one can reduce c/e while increase 
c3/V(c^e). 

Our analyses so far do not depend on the whether the inflation is caused by large 
field or small field, but when it comes to the quantitative analyses the differences are 
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interesting. For large field inflation e ~ r; while for small field inflation e <C ?/. For 
example, previous numerical works focus on the large field quadratic potential, and give 
the best fit jV ~ c ~ 0.2e ~ 0(10"^). Brane infiation is a small field infiation Q and 
e is much smaller, e ~ 0(10^^) or even 0(10"^^). Similarly the field range d (in units of 
Mpi) is also much smaller. So brane inflation is sensitive to very tiny steps present in the 
potential. This can be potentially used as a sensitive probe to the fundamental theory. 



5. Steps in IR DBI Inflation 



The DBI-CS action is Eq.(2.4) with V{<^^ = Vb — \m}<^ . The branes start from the tip of 



the throat and end at the UV end of the throat 4>cxiA = RbV^bTs = y/^/RB- After that, 
branes will quickly flnd antibranes and annihilate. 

In the relativistic (DBI) regime with t <^ —H~^, the attractor solution is 

This zero mode evolution relies on the background warped geometry, so it is reliable in 



regions where the back-reaction deformation from Hubble expansion is small [36, 38 1. This 
still leaves at least two interesting regions for primordial fluctuations: the field theoretic 
region (corresponding to smaller A'^e and shorter scales) and the stringy region (correspond- 
ing to bigger N^. and larger scales). The spectral index transitions from red to blue []36| , |3l| . 
Detailed estimates can be found in Ref. p5| . 

We now consider the effects of sharp features in the warped geometry on this IR DBI 
model. The duality cascade will leave sharp features on the warp factor T{(j)) but not the 
inflaton potential V{(p). Since the WMAP window spans a few e-folds, corresponding to 
branes rolling across Al « 3gsM/N^^^ duality cascades in a Klebanov-Strassler throat. 
In this paper we will mainly treat the steps as individual sharp features, namely they are 
well separated in terms of e-folds and sufficiently sharp to generate observational signature 
individually. For this purpose, we would like to consider the situation where 

gsM < Af . (5.2) 

In the multi-throat brane inflation scenario, the number of inflaton branes generated in B- 
throat is roughly determined by the flux number M. Since the number of inflaton branes 
is constrained to be around 10^ ~ 10^ |35], the condition ( ^.2| ) implies a small gs ^ 10~^. 



It is also very natural that gg takes value larger than 10 ^. In this case, the branes 
come across many steps within a single e-fold. We will discuss this case in Sec. ^.5| . 

5.1 The Properties of the Steps 

Here we like to make a crude estimate on the positions, fractional heights and widths of the 



steps, along the line of Sec. 4.1 but with more detailed numbers. The input is the position. 



fractional height and width of the step at / ~ 20 as well as the position of the assumed 
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step at / ~ 2. Let r = rg — Rb be the edge of the throat and r = at the bottom of the 
throat. They are related by the warp factor /i(</>f,) = hs, 

?'b~ro/is~roexp -- — — . (5.3) 



We are interested at regions where r;, <C r <C tq. Thus, given tq or r;,, we can determine 
ah the other transition values r^. Here p measures (starting at the bottom) the step in the 
duality cascade. In the large -fC S> p S> 1 limit, the correction is small, so we treat the 
steps as perturbations. The sharpness of the steps are controlled by the coefficient dp. A 
step becomes infinitely sharp as dp 0. We shall leave dp as free parameters. Note we 
can introduce similar spreads to the running of the dilaton. 

The position (pp of the steps can be computed. The steps are at (p = (j)p, where 

ln-^~ ip —. 5.4) 

Same as Eq.( |4.2| ) in the slow-roll case, we have 

-dNe ^ dink c::^ dlnl c::^ Hdt . (5.5) 

Suppose the pth step is at Ip (or (pp), because (p ~ —^/XB/t, we get 

(t>p tp-i iVo-ln/p_i 



6p_i tp iVo-ln/j 



(5.6) 



p 



where we have used the fact that H is approximately a constant, and A'^o is the number of 
e-folds (at the largest CMBR scale) to the end of DBI infiation, for example Nq ~ 38 |p5|| . 
On the other hand, duality cascade relates the position of the adjacent steps 

Suppose the feature at Ip^+i ~ 20 is due to such a step, and suppose the feature at Ip^ ~ 2 
is also due to a step (not just cosmic variance). That is, pi labels the initial or the first 
observable step. Then, taking A^o = 38, we see that the next 2 steps should be at 

lp^+2 ~ 170 lp^+3 ~ 1300 . (5.8) 

We also get a constraints on the microscopic parameters 

gsM ^ 33 . (5.9) 

Now we estimate the value of pi by looking at the dip in the CMB power spectrum at 
/ 20. The power spectrum is 

In IR DBI inflation, the Hubble scale is dominated by the constant y/Vo, so it does not 
change across a step in the warp factor. However, since (p closely tracks the speed limit 
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(f)^ ~ h'^ {(/)), it is most sensitive to the step in the warp factor. (/>^ increasing by a factor 
of (1 + 26) across a step decreases P-ji by a factor of (1 + 2b). The CMB data shows that 
around / ~ 20, there is a dip in the power spectrum by about 20% in fe-space, which gives 
26 ~ 0.2. Using we have 

^ * 0.2 . (5.11) 



Svr J {pi + 1)3 

Given that ggM ~ 33 to fit the spacing of the steps, we immediately get 

Pi = 2. (5.12) 

With input of 4 quantities : the position / ~ 20, the fractional height (size) AT/T ~ 0.2 
and the width A/p ~ 5 of the 2nd step as well as the position (at / ~ 2) of the first step, 
we can use A/ oc / and AT/T oc g^^^ complete set of predictions: 



V 


I 


AT/T 


Alp 


2 


~ 2 


~ 0.7 


~ 1 


3 


~ 20 


0.2 


~ 5 


4 


~ 170 


~ 0.08 


~ 40 


5 


~ 1300 


~ 0.04 


~ 260 



Note that AT/T ~ 0.7 and AT/T = 0.2 for the first two steps are probably too big to be 
treated as a perturbation. We shall take this to mean that the size of the first step can be 
very big. 

We should point out that the warp factor Eq. is a good approximation only when 
p ^ 1. This condition is also necessary to ignore the running of the dilaton. So p = 3 is 
already too small. Moreover, as we will see in Sec. the glitch in the power spectrum 
around / ^ 20 may be too large to be explained by a step feature in the warp factor due 
to the duality cascades. 

5.2 Qualitative Analyses around a Single Step 

We consider a sharp step in the warped geometry T(r). We parametrize the size of the 
step as 26 = AT/T and the width as A(j) = 2d. We use the tanh function interpolation 
between the two sides of the step, i.e. 



^4 



1 + 6 tanh 



(5.13) 



d 

5.2.1 The Evolution of the Sound Speed 

We use the case of positive 6 to illustrate the evolution of the sound speed, the formulae 
are the same for negative 6. The direction of the step is such that, for positive 6, T steps up 
as r increases, i.e. the warp factor h oc T^/^ (or the speed-limit h'^) increases as r increases. 

The first stage is when branes move across the step, where the speed-limit increases 
suddenly. Since the time that the branes spend to across the step is very short comparing 
to the Hubble time (see Appendix 1^) , during which we can ignore the Hubble expansion 
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and the acceleration from the potential. So across the step, the behavior of the branes 
is approximated by the Lagrangian L = —T{r)\J\ — (^P^ jTir). The corresponding energy 
T7 = const, where 7 = l/c^ = (1 — (jP' /T)~'^^'^ . Hence a sharp change in T results in a 
sharp change in c^, and both changes are small. 



Ac, 



AT 



2b . 



(5.14) 



A small change in Cg also means that the velocity of the branes will closely track the change 
of the speed-limit in step. 

The second stage is at the end of the step, where the previous change in Cg results in 
a deviation from the attractor solution. To see this, we note that in IR DBI inflation, the 
attractor solutions have 7 « /3-A''eDBi/3 and A'^eDBi ~ HR/h, so Cg « 2>h/ (3HR. A fractional 
change in the warp factor h = r/R means AR/R = —6/2, which implies that the attractor 
solutions has a fractional change 

Ac, 



b . 



(5.15) 



So the increase in Cs due to the first stage ( ^.14| ) over-shoots the new attractor solution, 
the branes will have to quickly relax to the new attractor solution. The Hubble expansion 
plays an important role here, so the energy Tj is no longer conserved. 



0.023 



0.022 



0.021 



0.020 - 





(t-fo)H 



Figure 3: Evolution of c^. Parameters are b — 0.1, step width AiVg = 0.05, /3 = 3, gsfn^ ^ = lO'^^, 
Nb = 10^ UB = 10^ nAh\ = 16. 

In summary, the behavior of Cg can be approximated by two step-functions side by 
side. In the first one, it jumps up from the original Cg to Cg + 2bcs, during which the 
branes quickly fall down the step in warp geometry. The time scale is determined by d 
and is typically much smaller than 1/H. The second one comes immediately afterwards, 
it jumps down from Cg + 2bcg to Cg + bcg, during which the branes quickly approach the 
new attractor solution. The time scale is of order 0{1/H). The width of the second step 
function is much larger than the first for narrower step (i.e. smaller d). Note that in both 
stages, the velocity of branes increases. This evolution is illustrated in Fig. ^. 
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Because the width of step is very small, such a small change of Cg happens in a short 
period. This gives rise to large s = Cs/{csH) and s/H. This turns out to be the primary 
sources to various observable signatures in power spectrum and non-Gaussianities. 

5.2.2 The Power Spectrum 

The power spectrum is determined by the mode equation Eq.(3.2). We now estimate 



various parameters in Eq. ( p.7[ ) in terms of the size, 6, and width, d, of the step feature in 
the warp factor. It is convenient to write d in terms of the e-fold that the brane spends 
moving across it, 

AiVe = HM ^Hi = —L= . (5.16) 

Note this is not the e-folds that the feature on CMB spans (denoted as Al previously) 
which typically includes the oscillations and spreads much wider. 

Using the continuity equation p = —3H{p + p), together with the energy density and 



pressure Eq. (|2.10D , we can solve for the evolution of s, 



2. . c,y T 



3(l-cf) + -^ + — (l-c) . (5.17) 



He, ' TH TH 



Here we can make more explicit the conditions under which the sharp features in warp 
factor dominate. On the right hand side of ( 5.17| ), to make the 3'^'^ term larger than the 



first, i.e. to ignore the spatial expansion, we require 3 < AT/{THAt) ~ 2b/ ANe] to ignore 
the 2^^*^ term, i.e. the sharp features in potential, we require CgV' < T'. Since the duality 
cascade leaves sharp steps on T{(p) but not V{(l)), this condition is easily satisfied. Under 
these conditions, the last terms in the above equation dominates, so we can approximate 

\/2i^~ O ( -r^ ) . (5.18) 



TH d \AN, 



This is consistent with ( |5.14| ). 

Denote X = ^'^ /2, so e = X/{M^^H^Cs) for DEI inflation. We get 



On the other hand, from the definition Cg = y^l — 2X/T{cj)), we get 

X f 2cls 
Yh ~TH~ H{l-cl) ' 

Plug (|5^ ) into (111), and use (|5.17D to get rid of f/TH, we have 



(5.20) 



fj = - 3(1 + Cs) - J!^^' + 2e 

Vr 

^CsS-{3 + j^) . (5.21) 
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Except for the first term, the r.h.s. of this expression is affected little, 0(b), by the feature. 
We use (5.21) to further estimate fj/H. Keeping the time derivatives on Cg, (p T, while 
ignoring those on V and H, we get 



(5.22) 



(5.23) 



We have used the fact that the brane always tracks the speed-limit, 

since <C 1 even at the sharp feature as we know from Sec. 5.2.1 . Eq. (|5.18|) is also used. 
We further estimate s/H, 



H V J T{(j)) \HdtH 



T" r'2\ 02 r if) 



Using Eq.(5.23) to eliminate the (p in Eq.(5.24), we get 

2c,e 



s 

H 



rpll 



ir'2\ fT" 



AN, 



(5.24) 



(5.25) 



where we have dropped the term T'^/T^, since T"/T » T'^/T^ as long as 6 < 1. 

We now can compare the amplitudes of various terms. Because both b and AA'^e are 
small, from ( |5.18| ) and ( 5.25 ), we have 



(5.26) 



From (15121] ) and ( |53^ ), 



So the only important contribution to z" / z is from s/H, and we have 

z \ 2H) \ T 

If we take T to be the form Eq.( |5.13| ), we can evaluate 

\ d J 

b 



z 



c.,e 



1 — 6-^5- sech^ 



1 



AiV2 



sech 



d 



tanh 



d 



(5.27) 



(5.28) 



(5.29) 



The feature in z" j z is dictated by the sech^ tanh term, and will affect the evolution of 
certain perturbation modes vy^. As moving across the step only generates ~ N^^^ jM 
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e-folds (See Appendix typically b/ AN^ ^> 1. The mode Vk will first see a dip in z" / z 
followed by a bump. Let us assume that the feature in z" /z shows up around conformal 
time Ts, then modes with c^fc^ ^ {z" / z)\r^ or c^fe^ <C {z" /z)\r^ will not be affected, because 
they are either oscillating well inside the sound horizon or have already crossed the sound 
horizon and got frozen. The major effect will be on the modes with c^k^ ~ {z"/z)\r^. The 
range of k affected by the sharp feature is determined by b/AN^. If Pfi starts seeing the 
feature at ko, the feature will disappear at 



Ak ~ koVb/ANe 



(5.30) 



5.2.3 Non-Gaussianities 

Although Ref. |41, 44 1 are only interested in non-Gaussianities without sharp features, 
the cubic expansion of the perturbations is exact and does not rely on the assumption that 
various inflationary parameters e, fj and s are small. For DBI inflation, among all the terms 
in Ref. |^T|, the leading term in the cubic action that is responsible for sharp features in 
non-Gaussianity is 



a^e d 



rj 



(5.31) 



In the absence of sharp features, a term such as (a^e/c^)CC^ contributes to the non- 
Gaussianity estimator Jnl ~ l/c^- Therefore as a rough estimate the term ( ^.31 ) con- 
tributes 
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The cases with large s are most interesting for non-Gaussianities. Using (5.21), 
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There is another term 
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(5.34) 
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in the cubic action of Ref. |41] that contributes 
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(5.36) 



which is also possibly observable. The term ( 5.34 ) dominates for the most interesting cases. 

The net observable effect will be a nearly-scale-invariant large non-Gaussianity of order 
0{l/c1) plus the scale-dependent (oscillatory) modulation of order 0{As/cs). Similar to 
the slow-roll case, we can write /l^l*"-^ ~ {Vb/cs){Vb/ANe). As we have shown, the 
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first factor determines the amplitude of the ghtch in power spectrum through b, while 



the second factor determines the extension of its osciUations through (5.30). So there is 
a simple relation between the effects of sharp feature on the power spectrum and non- 
Gaussianity. For a same single visible glitch in power spectrum, the oscillatory amplitude 
of the bispectrum in DBI inflation is larger than that in slow-roll inflation mostly due to 
the factor of 1/cs which can make f^'i'^^'^ larger than 0(10); while for glitches that have 
long oscillatory extension (which can also be made invisible in power spectrum with tiny 
b), the associated bispectrum can be very large with tiny ANg. The estimates of b and 
AA^'e in Appendix |^ and |^ suggest that the well-separated steps in the duality cascade 
tend to fall into the latter category. We will later see independent evidence from fltting 
the power spectrum glitch to WMAP CMBR data. 

5.3 An Analytical Approximation for a Single Sharp Step 

In this subsection, we give some analytical treatment of the effect of sharp feature on power 
spectrum. The full analyses rely on numerical calculations. Nonetheless we can see how 
several properties emerge and we summarize them in the end of this subsection. 



In Section |5.2.1| , we see that the sound speed undergoes two adjacent jumps. Here 
we simplify it by just considering the overall net effect, namely we approximate it by one 
step- function from Cs to Cg = bcg. We approximate the transition to be inflnitely sharp so 
that an analytical analysis is possible. As we showed, the width of the flrst step function 
can be made very small as we reduce d, while that of the second is much wider and of order 
0{1/H) model-independently. So this approximation is only qualitative. 

We denote tq to be the moment inflaton passes the feature, and look at the solutions 
of the equation of motion ( ^.2| ) at both sides of tq. At r = rj", 

Vk = vi\ - . (5.37) 

At r = T+, 

Vk = Civil+ + C2V2I+ . (5.38) 
Here we have denoted the two linearly independent solutions as 

_ mH ikcsT)e-''"'^^ , (5.39) 

.,.-^(l-^.c.r)e--. (5.40) 

For r < To we only have vi because the vacuum is dominated by the Bunch-Davis vacuum. 

These two solutions are matched according to the boundary conditions at tq. These 
are contributed dominantly by the sharp change in z" / z in the last term of ( |3.2D . From the 
previous subsection, we know that the change in z = a\/2e/cs is dominated by the change 
in Cg. As mentioned we approximate 

Cs = Cs\- + bCs\-6{T - To) , (5.41) 
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where 6 is the step- function. So 



z- (r) , r < To 

zq - bzo9{T - To) , r = To 
z+(r) , r > To 



(5.42) 



where zq = z_(ro). Around tq, we then have 
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, d ^, , Az' ^, , 
-b—5{T - To) + (5(r - To) + 

UT Zq 



(5.43) 
(5.44) 



where the dots denotes terms that do not contribute upon integration over the infinitely 



smah region around tq, and Az' 



Because z' = aHz{l + r]/2 — s), we have 



Vk, 



z'li = aHz\'^. Integrating ( ^.44 ) once and twice, we get the two boundary conditions for 



VkL+ - VkL- 

^0 ^0 
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and 
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Matching (|5.37|) and (|]3|) we get 
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where we denote xq = Tokcs\^-. The power spectrum is 
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where Pt^q = ^^^/(Svr^eCsl^-). For small A;, xq —>■ 0, so ^ P-jio as expected. For large k, 
P-jz is oscillating. The fact that such oscillation extends to infinitely large k is the artifact 
of approximating the change in Cs to be infinitely sharp. Such an approximation makes 
the potential barrier z" /z infinitely high, and introduces the V2 component to all k modes. 
Realistically, the relaxation takes finite time AA'e, and the potential barrier in z" /z is 
finite. For modes with high frequency kcs ^ z"/z, it will not see the barrier in z" /z and 
the V2 component is absent. Therefore, for the mode 



iz" , Vb 



(5.50) 



the oscillation is damped away. Here /cq = aH/cg is the mode number crossing the horizon 
at Tq. 
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Figure 4: The approximate power spectrum, P-ji, as a function of xq = kcgTo, where tq is the time 
of the feature, due to a function step with b = 0.1. 

This power spectrum is illustrated in Fig. ^. The fact that the dip appears first for 
positive b can be seen analytically by taking the bxQ <^ 1 limit, 

Pn^ Pno(^l- 2b sin 2x0^ . (5.51) 

Here we summarize several general properties from this analytical approximation. 
Firstly, the oscillation behavior in power spectrum can be understood as being due to the 
change of vacuum caused by the sharp feature. Such a change introduces a second (neg- 
ative energy) component of in addition to the original Bunch-Davis (positive energy) 
component. The superposition results in a oscillation in k-space. Secondly, the matching 
of the two set of solutions are done at the location of the sharp feature xq = kcsTQ = k/ko, 
where ko = {—tqCs)~^ is the wave number crossing the horizon at the moment of the 
feature tq. The phase 2xq gives the oscillation wavelength ~ ko/2 in power spectrum 
model-independently. This implies that, in log k coordinate, the oscillation behavior looks 
similar for sharp features at different k^. Lastly, for a feature with finite width, the high 
frequency modes in the vacuum change adiabatically and remain in the Bunch-Davis vac- 
uum. So the effect of the feature decay away for those modes. These conclusions apply for 
both the DBI and slow-roll case. 

5.4 Numerical Analyses and Data Fitting 

We have analyzed qualitatively the behavior of the power spectrum around the step. How- 
ever, due to the non-trivial behavior of z" /z, the mode equation is hard to solve analyt- 
ically in general. To better understand the effect of the sharp feature, we numerically 
solve the equation of motion Eq.(2.13), evolve the mode equation Eq.( |3.2D starting from 





k{/Mpc) 



Figure 5: In the IR DBI scenario, we sliow tlie power spectrum, P-jz, when there is a sharp step in 
the warp factor. For the same step width AiVe, we show three cases with different step size b. The 
amplitude of the first dip and bump increases as we increase the step height b. A bump appears 
first for 6 < while a dip appears first for 6 > 0. 



when the mode is deep inside the sound horizon (Cg/c^ ^ ^" / ^) until the mode is frozen 
when stretched out of the sound horizon (Cg/c^ <C z" /z), and evaluate the power spectrum 
numerically using Eq. (|3.4|) . The numerical calculation also has the advantage that it is not 



necessary to assume certain simplification conditions in Sec. 5.2, such as b/ANe ^ 1. Our 
major results are shown in Fig. |5| and Fig. ^. 

Fig. ^ shows the step feature in the power spectrum Pti for steps with different height 
6, but same width AiVg. We see clearly that in the power spectrum, a dip appears first 
with positive b and a bump appears first with negative 6, which agrees with our analytic 
result Eq.( 5.4g| ). However such a sequence gets less clear after projecting from /c-space to 
l-space. 

In Fig. ^ we explore the step feature with different step width AA'^e- We see that the 
oscillation amplitude in Pji is quite insensitive to the step width d. By the same argument 
in our slow-roll analysis, the steps are sharp so that AN^ <C 1. The Cg and Cg have reached 
their maximum values and so is the amplitude of the oscillation which is controlled by b as 
in (l5l9|) . 

In Fig. we perform a fit to the WMAP data. With one single step in the warp 
factor, we try to fit the / ~ 20 dip in WMAP temperature anisotropy. Here, instead 
of performing a full MCMC search for the best fit model, we only show an example to 
illustrate how the local step feature improves the quality of data fit. The example we 
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k{/Mpc) 

Figure 6: Same as Fig. ^, with tlie step size b fixed, we show two cases with different step width 
ANe- Since in IR DBf model, the typical step width corresponds to AN^ <C 1, the dip and bump in 
the power spectrum is not sensitive to d. The width of the step affects only the range of oscillation 
in P-ji. 



show in Fig. ^ has the IR DBI model parameters ub = 6761, Nb = 4.315 x 10^,n^/i^ = 
0.01035, mj^y"^ = 2.567 x lO^^Mp,/? = 4.021, and the step parameters are AN^ = 0.105 
and b = —0.35. The step model, with likelihood C, — 21n£ = 5347.16 when compared 
to WMAP 3-year temperature and polarization data @, ^, 45, 46], has a better fit to the 
data, Ax^ = -2.76, over the best-fit scenario with constant running of the spectral index, 
also shown. In light of the four extra degrees of freedom^ determining the cascade power 
spectrum in comparison to the fiducial running model, however, this improvement in is 
not statistically significant. 

We also find that a positive b is not favored by data, so the first dip in the fc-space 
does not necessarily transform into a clear dip in the /-space. Numerically we find that to 
fit the I ~ 20 dip, we need \b\ ~ 0.3. However, if we take b = +0.3, it gives too much power 
on large scale, more than that is allowed by data. 

The reason that the glitch appears to be less sharper than that in slow-roll case is 
the following. In DBI inflation, a step in warp factor not only causes oscillation in power 
spectrum, but also changes the asymptotic speed-limit of the inflaton. So in addition to 
glitches, it also introduces a step in the power spectrum. Both amplitudes are controlled by 
the relative height of the step b. It is more difficult to generate a sharp glitch while keeping 



'^Three of them specify the location, height and width of the step, another one is from IR DBI (versus 
LCDM). 
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Figure 7: The CMB temperature power spectrum (top) and initial power spectrum (bottom left) 
in the presence of a step in the warp factor (full black line) and a best-fit scenario with constant 
running in the spectral index (red dashed line) . The presence of large scale features in the cascade 
spectrum (bottom right) lead to an improved fit with the WMAP 3-year data of Ay^ — —2.76 over 
the model with constant running in the spectral index. 



the latter effect compatible with data. While in slow-roll case, asymptotic velocities of the 
inflaton is not affected by a local step in potential. In this sense, for DBI inflation the 
glitches in power spectrum can be made sharper if one introduces a bump (or anything 
that does not change the asymptotic behaviors of the warp factor) instead of a step in the 
warping. 

We should also point out that |6| ~ 0.3 gives a step much larger than that given by 
Eq.(^.g|) ^, which typically goes like with p ^ 1. This means that if we take the steps 
generated by the duality cascade, most probably we will not see any observable effects 
in the power spectrum. However, as we estimated in Eq. ( 5.34| ), the non-Gaussianities 



associated with the step can be very large depending on parameters in Eqs. ([A.5|) and 



*It is possible that duality cascade will generate large steps if we take p ~ 1. However, the perturbative 
analysis of the warp factor breaks down with p ~ 1, and we do not have any analytic control in that regime. 
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(IClol ). This can leave distinctive features detectable by experiments. 



5.5 Closely Spaced Steps 

We have been mainly concentrated on features that are sharp enough and well-separated 
(by at least one e-fold) . We also see that there are parameter space where the feature is not 
sharp enough to show up in the CMB as observable signatures. In this case, the step is no 
longer important individually; but if there are parameter space such that, within one e-fold, 
the branes can go over many such small steps, a different kind of large non-Gaussianities 
can be generated through the "resonance mechanism" ||4^. These closely spaced small 
steps induce small but high frequency oscillations of the background evolution in e, fj and 
s. This oscillation can resonant with the mode functions when the mode is still within 
the horizon. As a result of this resonance, the non-Gaussianity integration picks up a large 
contribution while the two-point function is only slightly affected. This non-Gaussianity 
has different signatures from the sharp feature case p^]. 

In slow-roll model, the condition to achieve this resonance mechanism relies on the 
shape of the potential. Naively, since we already have a series of steps, we can increase the 
tilt of the potential to make the branes roll faster so they come across many steps in one 
e-fold. Realistically, one needs to make sure that such a tilt will not significantly reduce 
the already-tuned total number of e-folds. 

In IR DBI model, as we have discussed at the beginning of Sec. |5|, there is actually a 
quite natural parameter space that such a condition can be satisfied, 1 > Qs > 10^^. We 
leave a more detailed study to a future publication. 

6. Remarks 

The analysis of the step feature in inflation in comparison with WMAP data has been 
carried out for large field models, specifically the quadratic chaotic inflationary model 
where the inflaton field starts at values much bigger than the Planck mass (0 > 15Mpi) 
and then decreases to zero. In brane inflation, the inflaton being the position of the brane, 
is bounded by the size of the flux compactification volume, and so typically cp <^ Mpi. 
This requires a generalized analysis, which is carried out here. In the IR DBI model, the 
behavior of the inflaton is further modified due to the DBI kinetic term. We have studied 
the properties and the signatures of the steps in warp factor and see that typically the 
effect on the power spectrum is much smaller, but non-Gaussianities can be much larger. 
So non-Gaussianities becomes a more important tool to probe the sharp features in DBI 
inflation. Another likely scenario is when the steps are too close to show up as sharp 
features in the power spectrum. Their presence is expected to show up in the bi- and 
tri-spectra. It is important to study this case further. 

We end with some remarks on two broader aspects on brane inflation. In this paper 
we considered brane inflation as an effective one-field model driven by a potential in the 
radial direction. The location of the brane in each of the compact directions (one radial 
and five angular) is represented by a scalar field - brane inflation is therefore multi-field 
in nature. Multi-field models can be decomposed into adiabatic and isocurvature fields, 
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respectively describing motion along, and perpendicular to, an evolutionary trajectory in 
field space, where sharp features, such as sudden turns, in the trajectory can give rise to 



the interconversion of adiabatic and isocurvature modes |48, 49, 50 1. The introduction of 
isocurvature perturbations can significantly modify the primordial power spectrum, and 
hence tight constraints on their contributions can be placed from observed CMB and large 



scale structure spectra |51, 52|. It will be interesting to see if these effects can be realized 



in brane inflation models. 

In brane inflation, there are a number of possible sources of non-Gaussianities: (1) 
Even for a single D3 brane case, the inflaton is actually a six component field, namely, the 
radial mode plus the five angular modes. The angular modes are massless in the simplest 
scenario, though in general they are expected to pick up small masses. This is a multi-field 
inflation, where non-Gaussianity may arise; (2) the DBI action; (3) cosmic strings and 
their cosmic evolution, and (4) step-like behavior. In this paper, we focus on the step-like 
behavior due to the warp geometry, and its possible detection. If more than one source 
of non-Gaussianity are present, it will be important to disentangle them, for example via 
bi-and tri-spectra. In this sense, brane inflation can be very rich. 
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A. The Warp factor 

Let us consider the structure in the warp metric h{r) that rescales masses, where 

ds'^ = h\r){-dt^ + a(t)2dx2) + h-\r){dr^ + r^dsj) 
For the KS throat, we have the approximate Klebanov-Tseytlin solution ||l5| 



h-\r) = ( N + '-^IHr/R) + 1/41 ) , (A.l) 



4r4 V 27r 
where, with the volume of T^'^ ~ S"^ x given by vTT^ = 167r3/27 



= Ana'^g.N/v = ^Jl^^ (A.2) 
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and M is the RR-flux wrapping while K = N/M is the orthogonal NS-NS-flux. The 
locations of duality transitions in the KS throat are given by 



Note that the constant piece N in ( [A.l| ) is determined by the boundary condition. We see 
that the effective Z)3-brane charge is given by 

N,s = N + ^^Hr/R) (A.4) 

so that at r = r^, A'eff = N — IM = K — I. The term with 1/4 factor is introduced to ensure 
that the warp factor /i^(r) is monotonic for R > r > r^- 

The size of the step at r = rp in the warp factor h{r) for large p < K (where r = r^+i 
is at the edge of the throat) is estimated perturbatively to be (with string coupling 
valued at the step), 



Ah{rp) h^{rp) - h^{rp) f 3gsM\ 1 



h{rp) h{rp) \ 27r J p^' 



(A.5) 



where hy{rp) is the warp factor at r > and h^{rp) that at r < rp. The widths of the 
steps are estimated in Appendix Q This series of steps leads to a cascading feature in the 
warp factor. So, in general, we expect this multi-step feature to be generic; as we approach 
the infrared (decreasing r), the step size grows (as p decreases). The spacings between 
steps are roughly equal as a function of In r for large p. 

In another throat, it is entirely possible that the leading Ah{rp)/h{rp) appears at a 
different order, say instead of It is also possible that Ahivp) /h{rp) can be negative, 
instead of positive, so h[r) is "saw-like" instead of a cascade. 

Here we write down the warp factor with steps given in |19] for a flow including only 



s number of steps starting from some location of duality transition at r = ri^] i.e, in the 
range r^^.i < r < Hq+s- The warp factor in this case can be written as 



h-\r)^h-\lo,r) + '-^ ^ [gsPiM 



l=lo 



Sir 



X l + tanh[^ ] , (A.6) 

2pf V qri 



where pi = K — (l — 1) and h ^(/qi ^) is given by 



277ra 



'•2 r Q„2 n/r2 



3giM' 



^-^(^o,r):.^^,.Ar + ^[ln(-) + -] 



^»-*^'|ln(-) + 2(l„(^))^ + ll)^-^l. (A.7) 
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Note that the size of the step in h ^(r) at r = r; is 

where A'cff = N — (l — 1)M = piM is the effective D3-branes charge for the flow from the 
{I — l)*'^ to the duahty transition point. 

B. The z"/z in Slow Roll Case 

In the presence of a feature hke the step in the potential, the motion of the inflaton is 
divided into three phases: attractor before the step, acceleration within the step, and 
relaxation after the step. During the attractor phase, the inflaton moves along with the 
slow-roll attractor determined by the flat potential. In the acceleration phase, the inflaton 
gets accelerated by the local steep potential due to the step, and deviate away from the 
attractor solution. In the relaxation phase, the potential is flat again and the inflaton 
motion relaxes back to the slow-roll attractor due to Hubble friction. 

Since Hubble friction is responsible for damping the deviation from slow-roll attractor, 
it is natural to expect that the relaxation time for cj) is of order or 0(1) e-fold, and the 
same for e, fj. However, the sharp change in z" /z takes place within a period much shorter 
than 0(1) e-folds. To see this, we first express z" /z in terms of inflation parameters. 



2a'H'^l + lfj + e) , (B.l) 



where in the last line we have ignored the e term since e ^ 1 in small field brane inflation 
models. The definitions of e and fj follow the main text and is defined the following way 



.2 _ 1 



We can use the following equation of motion to eliminate the term d^(p/dt^ in 

dV 

cP + 3H(I) + — = . (B.3) 
aq) 

We then get 

.2 3 V'U) 3 

«^ = -^ + 3.--,.-^-j,. (B.4) 
Combining ( |B.1| ) and (B^), we get 



The result ( [B.5D is very useful in determining the relaxation time of z" / z. We see that the 
sharp feature in z" / z vanishes once V"{(f)) goes back to normal, which only requires that 
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(j) has moved across the step. Since the potential is very steep at the step, it only takes 
AA'^e <C 1 e-folds for (p to move across. So z" / z relaxes much faster than e, 77 and When 
z" /z goes back to the attr actor behavior, e, rj and still needs a few more e-folds to relax. 
Both are important to observational quantities. 

C. Estimate the width of the steps 
C.l The width from Seiberg duality 

The running of the two coupling constants gi and 52 behave as 

The constant go and the scale A are chosen such that at /j, = A, gi = g2 ^ go- In terms of 
the string coupling gs, go = ^ngse^, because 

87r2 87r2 _ 27r 
9i 92 5s e* 

The duality cascade happens when one of coupling, for example gi, becomes large. We 
denote such a coupling to be gsD, which we expect to be of order one. The corresponding 
scale is 

liSD = Ae^^^~^^ . (C.3) 
The scale where one of the coupling hits infinity is 

/iloo = Ae^ , (C.4) 

1x200 = Ae ^ . (C.5) 

So jioo is the location of the transition wall, and jisD characterize its thickness. Therefore 
the width of the step is 

^^l■SD = 2|//ioo - MSdI , (C.6) 
and the separation between the two steps is 

Imoo - l^2oo\ ■ (C.7) 

Their ratio is 

A o_2 



1 r « ^— 5se • (C-8) 

l/^loo — ^2oo| gsD 

In terms of the distance r, the separation between the step I and / + 1 is 2nri/3gsM. So 
the width of the step is 

IGvr^ r; 
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In IR DBI inflation it is also convenient to write such a width in terms of the e-folds 
that brane spends crossing the width of the step. From the attractor solution r ~ B? /\t\ 
and the fact that the speed-limit does not change much across the step, we get 



AiVf^-S^^, (CIO) 



IGvr^ AT, 

where Nf^ is the e-fold to the end of DBI inflation at the location of the step. 



C.2 The width from multiple brane spreading 



In the multi-throat brane inflation scenario |13|, inflaton branes are almost always generated 



in a large number. A typically generation mechanism can be the flux-antibrane annihilation. 
Once they are created at the end of the annihilation process, they can have different 
velocities. This introduces a spread in the mobile brane position. In this subsection, we 
estimate such a width. We express it in terms of e-folds 

The time-delay between the fastest brane and the static brane is the time period 
that the potential accelerates the brane from cf) = Q io the speed-limit (jP' /^/N. Such an 
acceleration can be approximately described by 

;f) = V'{^) =m^cp . (C.ll) 

Consider initially, at t = 0, 4>min = and 4>min = VTarmin = VNhmin/R- At t = At, 
(j) = (j)^/\fN. Here N and R are the charge and scale of the throat, respectively. We can 
therefore solve for At assuming mAt <C 1. Further, using m ~ ff, we get 

HR 

For a long throat, the minimum warp factor is typically determined by the Hubble defor- 



mation as a result of the dS space back-reaction. In Ref. |36, 38], it is estimated to be 

hmin ~ HR/ y/N. So 



Comparing ( |Clo| ) and (|Cl^) , we see that generically AiVf-^ > ANT''^'^ , although the 
relation can change for special parameters. The common feature is that both of them are 
almost invariant for different steps. 
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